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1 . INTRODUCTION 


Recently,  a number  of  papers  have  discussed  ferromagnetic  and  anti- 
ferromagnetic phases  of  Ionic  solids  containing  a rare  earth  element  as 
one  of  the  major  constituents.^"**  In  one  of  these  papers,^  a second 
rare  earth  Is  Intentionally  substituted  In  a small  percentage  to  be  used 
as  a probe  to  Investigate  the  Internal  fields  In  these  solids  at  temper- 
atures below  the  phase  transition.  The  data  were  analyzed  on  the 
assumption  of  an  effective  magnetic  field  Inside  the  sample.  It  was 
found  that  no  single  value  of  the  effective  field  could  be  chosen  to  fit 
the  Zeeman  splitting  observed  on  the  substituted  Impurity  Ion. 

In  some  of  the  papers,  the  experimental  data  were  Interpreted  In 
terms  of  an  effective  magnetic  field  derived  from  a sum  of  the  magnetic 
dipoles  taken  over  the  lattice.  We  show  that  these  terms  are  only  the 
lowest  order  of  a particular  multipolar  expansion  caused  by  the  Inclu- 
sion of  the  finite  extent  of  the  4f  electrons  of  the  rare  earth  Ions. 

In  this  report,  we  first  derive  a convenient  form  for  the 
Interaction  Hamiltonian  of  an  orbital  electron  In  the  presence  of  a 
remote  dipole  and  then  cast  this  Interaction  Into  Irreducible  tensor 
form.  The  single-electron  matrix  elements  are  calculated,  and  the 
results  are  used  to  cast  the  Interaction  Into  unit  spherical  tensor  form 
for  the  N electron  configuration  (nt)*^.  The  coefficients  of  the  unit 
spherical  tensors  are  shown  to  be  particular  lattice  sums  Intimately 
related  to  those  obtained  In  the  point  charge  model  of  Stark  splitting 
of  Ions  in  a crystal. 


2.  FUNDAMENTAL  INTERACTION 
2.1  Magnetic  Dipole  Field 

We  are  concerned  here  with  the  effect  of  an  Inhomogeneous  mag- 
netic field  due  to  external  dipoles  on  the  energy  levels  of  a paramag- 
netic Ion.  Thus,  we  are  led  to  consider  the  fundamental  Interaction 
between  an  electron  and  an  Inhomogeneous  magnetic  field.  In  this 
report,  the  sources  of  the  Inhomogeneous  field  may  be  considered  pure 
dipoles. 

^Donald  Jean  Randazzo,  PhD  Thesis,  ’ Johns  Hopkins  University, 
Baltimore,  MD  (1966)  (University  Microfilms,  Ann  Arbor,  MI, 
No.  66-12,514). 

2p.  J.  Becker,  M.  J.  M.  Leask,  and  R.  N.  Tyte,  J.  Phys.  C,  5 (1972), 
2027.  ~ 

3j.  C.  Wright  and  H.  W.  Moos,  Phys.  Lett.,  29A  (1969),  495. 

**H.  G.  Kahle,  L.  Klein,  G.  Muller-Voght , and  H.  C.  Schopper,  Phys. 
Status  Solidi  (b) , 44  (1971),  619. 


The  magnetic  field  at  the  point  f due  to  a dipole  of  moment  M 
situated  at  the  point  R is  given  byS 


. 3[M  • (R  - r)1  (R  - r)  M 

|R-r|'  ‘|R-r|3 


(1) 


It  is  well  known  that  the  magnetic  vector  potential  corresponding  to 
such  a field  is  given  by® 


I = . . (2) 

This  vector  potential  may  be  written  in  the  form* 

A = - [M  X grad  ^(r)],  (3) 

where 

i){r)  = ^ . (4) 

|R  - r| 

Equation  (4)  may  be  expanded  in  spherical  tensors: 

" k^q  ^kq  (5) 

where 

(6) 

and  where  Ykq  is  the  ordinary  spherical  harmonic.® 


®iV.  K.  H,  Panofaky  and  M,  Phillips,  Classical  Electricity  and 
Magnetism,  Addison  Wesley,  Reading,  MA  (1955),  127-128. 

E.  Rose,  Elementary  Theory  of  Angular  Momentum,  John  Wiley  and 
Sons,  Inc.,  Wee  York  (1957). 

*Here  and  in  what  follows,  we  use  the  notations  "grad,"  "curl,"  etc., 
to  represent  an  operator  acting  on  only  the  function  directly  on  the 
right  of  the  operator,  whereas  the  operator  "V"  acts  on  everything  to 
the  right.  For  example, 

V((i  = (grad  (j.)  + ((i7. 
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Here  the  problem  1s  to  determine  the  Interaction  between  a 
paramagnetic  Ion  and  a lattice  of  dipoles.  Thus,  we  must  sum  the  fields 
of  equation  (1)  over  all  dipoles  In  the  lattice.  Mere,  we  make  our  most 
crucial  assumption:  that  the  moments  of  all  the  paramagnetic  Ions  In 
the  lattice  are  equal;  we  allow  for  the  possibility  that  the  signs  might 
vary,  to  allow  for  an  antiferromagnetic  lattice.  With  this  assumption, 
we  may  use  equation  (3)  for  the  vector  potential,  where  *(?)  Is  now 
given  by 


1>(r)  = I 
k,q 


'kq  ^kq 


(r)  r' 


where 


N 


kq  ■ I 

^ paramagnetic 


C (R,) 


ions 


(R,.) 


WT 


(7) 


(8) 


kP1 


and  where  the  factor  (-)  Is  considered  to  take  account  of  the  varia- 
tion In  sign  of  the  dipoles.  The  quantities  Nkq  are  called  the  magnetic 
crystal  field  parameters,  by  analog  with  similar  quantities  which  occur 
In  considering  the  Interaction  of  Ions  with  the  Inhomogeneous  electric 
field  of  crystal  lattices.^ 

2.2  Correct  Expression  for  Interaction 

The  expression  for  the  Interaction  between  an  electron  and  a 
magnetic  field  Is® 


Hm  = 1;;:  (P  • A + A . p)  + (?  • curl  A) 


(9) 


If  Into  equation  (9)  we  substitute  equation  (3)  for  the.^ctor  poten- 
tial, tak1 
find  that 


_^Ct( 

tial,  taking  account  of  the  fact  that  the  divergence  of  A Is  zero,  we 


Interactions 

Laboratories 


^N.  Karayianis  and  C.  A.  Morrison,  Rare  Barth  Ion-Host 
1.  Point  Charge  Lattice  Sum  in  Scheelites',  Barry  Diamond 
TR-1648  (October  1973). 

®L.  J.  Schiff,  Ouantum  Mechanics,  3rd  ed.,  McGratf-Hill  Book  Co.,  New 
York  (1968),  478. 
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= 2yg  [1M  • (grad  (|»(r)  x 7)  + (M  • grad)  (s  • grad  <t)(r))]  . (10) 

In  equation  (10),  yB  * ef)/2mc,  and  M Is  the  dipole  moment  of  an  external 
dipole.  The  quantity  t Is  the  spin  angular  momentum  operator  of  the 
electron. 

We  now  cast  equation  (10)  Into  Irreducible  tensor  form.® 


3.  IRREDUCIBLE  TENSOR  FORM  OF 

We  shall  not  consider  the  full  Interaction  Hamiltonian  of 
equation  (10),  but  rather  we  shall  break  It  Into  portions  according  to 
the  decomposition  of  equation  (7).  Thus,  equation  (10)  becomes 


H„  - Nj,  H,,  . (1,) 

where 


Hj^q  = 2yg  [IM  • (grad  (t)|^q  x v)  + (M  • grad)(s  • grad  (j)|^q)]  (12) 


and  where 


(r) 


(13) 


We  now  cast  equation  (12)  Into  Irreducible  tensor  form.  Using 
equation  (A-26)  (app  A),  we  find  that  the  second  term  of  equation  (12) 
becomes 

Hk  (spin)  = 2yg  I M s_^(grad)  (grad)^  <^k 

n.v  ^ 

= 2ug  [k(k  - 1)  (2k  - 1)  (2k  - 3)]*^  r*^’^ 

X <k-2(q+y+v)l(-y)  I k-l(q+v)><k-l(q+v)l(-v)  I k(q)> 

P.v 

^k-2,q+y+v  ^-v  • (14) 

®W.  E.  Rose,  Elementary  Theory  of  Angular  Momentum,  John  Wiley  and 
Sons,  Inc.,  New  York  (1957). 
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It  1s  convenient  to  recast  equation  (14)  Into  a form  which  reflects 
explicitly  that  It  Is  a tensor  of  rank  (k-1)  In  the  total  angular  momen- 
tum space  of  the  electron.  It  may  be  recast  by  recoupling  Its  Clebsch- 
Gordan  (C-G)  coefficients  In  equation  (14)  according  to  equation  (A-13) 
(app  A): 

<k-2  (q+vi+v)l{-u)  I k-l(q+v)><k-l(q+v)l(-v)  | k(q)> 

= (2k-l)  W(1  k-2  k 1;  k-1  k-1) 

X <k-2(q+u+v)l(-v) 1 k-l(q+u)><k-l(q+vi)l(-w) 1 k(q)>  , (15) 


But 

W(1  k-2  k 1;  k - 1 k - 1)  = ^k^  1 


Therefore,  we  may  rewrite  equation  (14)  as 

Hkq(spln)  = 


= 2ug  [k(k  - 1)  (2k  - 1)  (2k  - 3)]*®  r'^"^ 

^ I <k-2{q+y+v)l(-v) I k-l(q+y)><k-l(q+y)l{-u) (k(q)> 


X C^,  , , s M 

K-2,q+u+v  -V  -y 


(17) 


Now  the  first  term  In  equation  (12)  Is  considerably  more  difficult 
to  reduce.  We  may  write 

Hkq(orbit)  = 2^/^  Pg  I (.)^  M <l(v)l(y-v) | l(y )> 

^ y,v  ^ 


X (grad  (j>,  ) 7 

^kq'v  y-v 


■i/2 


= 2"'  y.  I (-)^  M <l(v)l(y-v)|l(y)>(grad  <D.  ) 


y ,v 


kq'v 


A 

9 


(18) 


Using  the  gradient  formula  equation  (A-26)  (app  A)  in  equation  (18), 
we  find,  using  equation  (13)  for  (grad  <(>|^p)v. 


3/2 


H^q(orbit)  = 2 (2k+l)<l(0)  klO)  | k-l(0)>pg  I (-I)'' 


y.v 


X <l(v)  l(p-v)  lUy)>  <l(v)  k(q)  lk-l(q+v}> 


* 'k-l,q+v  If  - 'H2kn)<l(0)kl0)|k-1(0)>  , 


B 


X I (-1)^  M_  <Uv)Up-v)  I Up)><Uv)k(qj  I k-l(q+v)> 


y,v,X 


Thus,  the  orbit  part  of  the  interaction  divides  into  two  parts;  we  shall 
consider  each  part  separately. 


The  first  term  in  the  interaction  is 


(2)  _ ,3/2 


% - 2 yg  (2k  + l)<l(0)k(0)lk-l(0)>  J 


y ,v 


X <l(v)l{y-v)  l(y)><l(v)k(q) I k-l(q+v)>  r*^"^ 


X C 


k-l,q+v  i,y-v  9r 


(20) 


.3/2 


2 yg  (2k  + 1)<1(0)  k(0)|k-l(0)> 


x^  <k-l(0)l(0)lk'(0)>  I (-)^^  M 


y.v 


X <l(v)l(y-v)|l(y)><l(v)k(q)|k-l(q+v)> 


X <k-l(q+v)l(y-v)|k'(q+y)>  r*^"^  |p 


(21) 
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By  recoupling  the  last  two  coefficients  and  summing  over  v,  we  obtain 


= 2pg  [6  (2k  - 1,-^  (2k  + 1)<1(0)  k(0)|k-l(0)> 

X <k-l(0)l(0)|k(0)>W(kl  kl;  lk-l)r'^‘^^<k(q+u)l(-p) | k(q)> 

u 

X C.  . M • (22) 

k,q+u  -y  3r 

Now  when  we  evaluate  the  C-G  and  Racah  coefficients  of  equation  (22),  we 
find 

HkJ,"^  = 2yg  [k(k+l)]’5  r*"-^  I <k(q+y)  l(-u)|k(q)>  M_^  . (23) 

We  now  come  to  the  third  term  In  the  Hamiltonian,  the  second  term  In 
equation  (19).  This  may  be  written  as 

H = - 4^0  (2k  + l)r'^‘^<l(0)k(0)lk-l(0)>  I I 

M o y,v,x 

X <l(y)k(q)|k-l(q+v)><l(A)l(y-v-x)ll(y-v)> ^IX^y-v-X 
=-4yg(2k+l)<l(0)k(0)lk-l(0)>r'^"^  ^ <k-l(0)l(0) | k' (0)> 

^ <1( v-X)l(y-v+x) |l(y)><l(v-X)k(q)| k-l(q+v-X)> 

* <l(x)l(y-v) I l(y-v+x)><k-l(q+v-x)l(x) | k' (q+v)> 

xC..„.  ^ (-)‘'M  . (24) 

k Vl+v  y-v  -y 

We  may  recouple  the  first  two  C-G  coefficients  and  then  the  last  two 
and,  finally,  sum  over  x to  obtain 

= -12y„[k(2k-l)]V-2  Y [(2k'+l)(2k"+l)]‘' 

X W(1  1 k-1  k'ak")W(l  1 k-1  k';  Ik")  <k"(q+y)l(-y) | k(q)> 

X <k'(q+v)l(y-v)|k"(q+y)><k-l(0)l(0)|k'(0)> 

xC.,  . t M . (25) 

K,q+v  y-v  -y  ' ' 
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We  now  repeat  the  expression  for  in  irreducible  tensor  form: 

”kq  ^kq  ^kq  ^ ^kq  ’ 

where 


Hkq'^  = 2yg[k(k-l)(2k-l)(2k-3)r*  I <k-2(q+u+v)l(-v)  | k-l(q+y)> 

y,v 


1^2  ..k-2 


x<k-l{q+y)l(-y)|k(q)>  C,^.2.q+y+v  ^y  ’ 


^kq"^  = 2yg[k(k+l)P  r"-^  I <k(q+y)l(-y)|k(q)>  M_^  ^ 


>1 


= - 12yg[(k(2k-l)f^  r''-^  { _/-)'''‘'''"[(2k'+l)(2k"+l)f* 

k j k 

^ W(1  1 k-lk;  lk")W(l  1 k-lk';  lk‘')<k-l(0)l(0) ( k ' (0)> 


X [ <k'(q+v)l(y-v)(k"(q+y}><k"(q+y)l(-y))k(q)> 

y ,v 

X M C.  , , £ 

-y  k ,q+v  y-v 


4.  SINGLE  ELECTRON  MATRIX  ELEMENTS 

We  consider  matrix  elements  of  the  operators  of  equations  (27)  to 
(29)  between  wave  functions  of  the  form 


) x('"c>  * * 


s'  ''""s'  ”m  m' 
s s 


First  of  all,  we  note  that,  1f  n=n'  and  t=t’,  the  matrix  elements  of 
equation  (28)  may  be  Integrated  by  parts  over  the  radial  variable  to 
obtain 


X Cl,  . M 

k,q+u  -VI 


(32) 


Thus,  according  to  the  Wigner-Eckart  theorem,  we  obtain  for  the  matrix 
elements  of 

<£m^'s  smg>  = -pg  (k  + l)[k(k  + 1)]^  <r'^“^>  6^ 

<-^llC,^||'C>  I M_^<k(q+p)l{-u)lk(q)><£(m^)k(q  + y)l£(m^'): 
u 

where 

nl^ 


(33) 


dr  r'^lR„^(r)p. 

0 


(34) 


The  quantity  <)i|(C|^[(ji>  Is  the  reduced  matrix  element  of  C.  (app  B). 


The  matrix  elements  of  H 


(3) 

kq 


are 


<£m^'s  IT’s' m^>  = -12Mg[k(2k-l)]^<r*^"^ 


>6m^m' 
s s 


I (-)'^'‘^'^"[£(£+l)(2k'+l)(2k"+l)f^<£|  1C  '1  l£><k-l(0)l(0)|k'(0)> 


k'  ,k" 


X W(1  1 k-lk  ;lk")  W(1  1 k-lk';  Ik")  [ <k' (q+v)l(y-v) | k"(q+p)> 

y ,v 

X •<k"(q+y)l(-y) lk(q)>  <£(m^)l(y-y) |£(m^+y-v)> 

X <£(m^+y-v)  k' (q+v) |£(m£' )>  . (35) 
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We  may  recouple®  the  last  two  C-G  coefficients  of  equation  (34)  and  sum 
on  V to  obtain 

<£m^'s  nis|H^J|^^|£m^s  m^>  = 12ijg[k(2k-l)£(£+l)(2£+l) 

m.  I (2k"+l)/(2k'+l)<k-l(0)l(0)lk'(0)>  . 

k',k" 

X W(£Uk';£k")W(l  1 k-lk  ;lk")W(l  1 k-lk';  lk")<£|  |C,^' | U> 

xl  M <£(m^)k"{q+ii)|£(m^')><k"(q+M)l(-p)|k(q)>.  (36) 

u 

To  determine  the  behavior  of  equation  (35),  we  must  examine  each 
term  separately.  The  quantity  of  Interest  In  our  analysis  Is 

t(k,k'.k")  = (2k"+l)[k(2k-l)(2k'+l)£(£+l)(2£+l)]^<k-l(0)l(0)|k'(0)> 

X W(£Uk';£k")W(l  1 k-lk;  lk")W(l  1 k-lk' ;lk")<£|  |C,^|  |^>.  (37) 

.By  evaluating  explicitly  the  quantities  occurring  In  equation  (37) 
according  to  equation  (A-17)  (app  A), 

t(k.  k.  k)  = y2  (k  + 1)  /k(k  + 1)  <t\ |C^| \t>  , (38) 

t(k,  k.  k-1)  = - jkg-g  - .ll(2-e  + -k  + 1)  <^||C,^||^>  , (39) 

t(k,  k - 2,  k - 1)  = - ^ ~ l^kl 


where  we  have  used  the  equation. 


V 


Therefore,  the  total  of  equation  (36)  is 


stn5>  = yglk+D/kTk+TT  <r'''^>6^  <£l|C^||£> 


s s 


X J M_^  <«,(m^)k(q+y)l£(nip><k(q+y)l(-y)  |k(q)> 


-ViB[k(2k-l)(2£-k+l)(2£+k+l)]'^  m'  '='^llCkll^> 


s s 


* 1 M_^<£{m^)k-l(q+y) |£(m^' )><k-l(q+y)l{-y) I k(q)>. 


(42) 


The  first  ter-m  of  equation  (42)  cancels  equation  (33)  exactly,  and 
therefore  we  have 

<£m£'s  mg>  = - yg[k(2k-l)(2£-k+l)(2^+k+l)]^= 


X M_^<£(m^)k-l(q+y)K(m^')> 


s s 

X <k-l{q+y)l(.y)|k(q)>  . (43) 

Finally,  for  the  spin  part,  we  have 

^£in^'s  ms'l^kq^^l'^£^V  " 2yg[k(k-l)(2k-l)(2k-3)s(s+l)]^<£|  |C,^.2l  |-e> 

X <r^  I M_  <k-2(q+y+v)l(-v) |k-l(q+y)>  <k-l(q+y)l(-y) jk(q)> 
y ,v 

X <£(m^)k-2(q+y+v’)|£(m^')><s(tng)l(-v)|s(m^)>  . (44) 


5.  MATRIX  ELEMENTS  FOR  EQUIVALENT  ELECTRONS 

We  now  consider  the  problem  of  N electrons,  each  with  orbital 
angular  momentum  t.  We  introduce  the  unit  tensors®  uj^  and  v^^,  which 
are  defined  such  that 


9g.  Racah,  Phys.  Rev,,  ^ (1942),  438, 
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<£'s' ||uN|£s>  = 6^’  655’  (45) 

and 

<^'s'||v^^|£s>  = 6^'  655'  . (46) 

k Ik 

The  quantity  u Is  a tensor  of  rank  k In  orbital  space;  v Is  a double 

tensor  of  rank  k In  orbital  space  and  rank  1 In  spin  space.  We  may 
write  the  Hamiltonian  equations  (27)  to  (29)  In  a form  using  the  unit 
tensors : 

HkJ'^  = 2uB[k(2k-l)(2k-3)(k-l)s(s+l)]’^<£||C^_2ll-^><r''‘^> 

I M ^^l^^2(q+y+v)l(-v)|k-l(q+y)> 

VI, V 

X <k-l(q+w)l(-y)|k(q)>.  (47) 


HkJ^’  = -UB[k(2k-l)(2^-k+l)(2^+k+l)]^^||C,^ll^><r''•^ 


X I M_^  Uq;J<k-l(q+vi)l(-y)|k(q)>  , 


(48) 


with 


H a U ( ^ ) + u (2') 

”kq  ^^kq  ^ ^^kq  ' 


(49) 


When  we  qo  over  to  N equivalent  electrons,  we  sun  the  one-electron 
Hamiltonian  equation  (49)  over  all  N electrons.  We  define  the  unit 
tensors  for  N electrons  as  the  sum  over  all  the  one-electron  tensors: 


U 


k 

q 


N k 

I 

1=1  ^ 


wlk 

Aq 


N 


(i). 


(50) 

(51) 
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Therefore  the  Hamiltonian  for  N equivalent  electrons  becomes 

"kqCl  ' "kV’  * C’  ■ 


(52) 


where 


= 2yg  [lc(k-l)(2k-l)(2k-3)s(s+l)]*®<£llCk_2lU><r‘''^> 

X I M-J-v'^q+y+v  <k-2(q+y+v)l(-v)  | k-l(q+y)> 
y.v 

X <k-l(q+y)l(-y) 1 k(q)>  (53) 


and  where 

Hkq"'^  = -yg  [k(2k-l)(2^-k+l)(2£+k+l)]^£llCklk><r'''^> 

X I M_  <k-l(q+y)l(-y)|k(q)>.  (54) 

When  matrix  elements  are  taken  of  equations  (54)  and  (53),  they  are 
expressed  in  terms  of  the  C-G  coefficients  and  reduced  matrix  elements 
of  the  unit  tensors.  In  particular,  in  a basis  |JMLS>,  where  L is  the 
total  orbital  angular  momentum,  S is  the  total  spin  angular  momentum, 
and  J = t + S,  we  have 

<ct'J*M‘L'S‘lHkJ^^laJMLS>  = 2yg[k(k-l)(2k-l)(2k-3) s(s+l)]^<£l 10^.2! 

X <r'^'^<a'L’S'|lV^’'''^||aLS>[(2s+l)(2S'+l)(2L'+l)(2k-l)]^ 

X X(SLJ;  S'L'J';  1 k-2  k-1)^  <k-l(q+y)l(-y) | k(q)> 

V* 

X <J(M)k-l(q+y)|J'(M')>  M_^  (55) 

for  the  spin  part  of  the  interaction  and 
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<a'J'M'L'S' |H,^(2')!aSMLS>  = -yg[k(2k-l)(2^-k+l)(2^+k+l)]^  <£||C,^||£> 
X <r'^’^><a'L'S'|  |aLS>[(2L'+l)(2J+l)]*®  6^^,  W(k-1  L'J  S;  L J') 


X I <k-l(q+y)(-y)|k(q)><J(M)k-l{q+y)|J'(M')>M_^  . (56) 
y 

The  quantities  a and  a'  In  equations  (55)  and  (56)  are  additional  quan- 
tum numbers  necessary  to  uniquely  Identify  the  states.  The  quantity 
X(abc;def;gh1)  is  the  X-coeffIclent  and  Is  Identical  to  a 9-J  symbol. 
The  manipulations  necessary  for  the  results  In  equations  (55)  and  (56) 
are  given  In  appendix  B.  The  quantities  <a'L'S' | | |oLS>  and 
<o'L'S' I ||oLS>  are  reduced  matrix  elements  of  the  unit  tensors  In 
the  Russell -Saunders  basis  and  are  tabulated  extensively.^^ 


6.  OPERATOR  EQUIVALENT  FOR  LOWEST  TERM 


We  now  consider  the  k = 2 part  of  the  Hamiltonian  equations  (47)  and 
(48).  Substituting  the  explicit  form 


<£||C2l|£>  = <^(0)2(0)|£(0)>  = - 

the  Hamiltonian  becomes 


lU+1) 
(2£+3)(2£-l) 


(57) 


^2q  = I'B  ^ ^-y  <l(q+y)l(-y)  |2(q)> 


(58) 


lOc.  w.  Nielson  and  George  F.  Foster,  Spectroscopic  Coefficients  for 
the  p”,  and  Configurations,  The  MIT  Press,  Cambridge,  MA  (1963), 
The  relations  between  the  reduced  matrix  elements  used  in  that  book  and 
in  our  report  are 

<a'L'5'|  |/|  |ar.S>  (o'L'S*||t/'l|oIs) 

and 


<a'L'S'| |aLS> 


2 r 21  + 1 

7J  l(2f  + 1)(2S’  + 1) 


lo'L'5' 
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We  now  prove  that  this  equation  Is  Identical  to  the  usual 
Hamiltonian  for  a constant  magnetic  field: 


H = ug  (£  + 2s)  • H . 
We  consider  a single  term, 

d).  = C,  r*^ 

^kq  kq  • 


(59) 

(60) 


and  calculate  the  field  according  to  equation  (3): 


H = - [7x(i^xV(t))]  = 2 'l  <l(v)l(y-v)  I l(y)><l(X)Uy-v-X)  1 l(p-v)> 

I'  ^ vX 

X M V 7 .C,„r2=30M„.  <1(0)2(0) | 1(0)><1(0)1(0) 10(0)> 

X V y-v-X  tq  M'^y 

X I <l(v)l(y-v)|l(y)><l(q+y)l(-q-v)ll(y-v)><l(-q-v)2(q)|l(-v)> 

V 

x <l(v)l(-v)l0(0)>  . (61) 

If  we  recouple  the  first  two  C-G  coefficients  and  explicitly  evaluate 
them  with  j = 0,  we  arrive  at  the  result 

H = (-)^‘  /e  M„ . <l(q+y)l(-y)|2(q)>  • 

y q+y 

This  result,  when  substituted  in  equation  (59),  gives 

H = ye  ^ ^^-y  • 


(62) 


(63) 


which  was  to  be  proved. 


7.  DISCUSSION  OF  RESULTS 

In  view  of  the  results  of  section  6,  It  Is  convenient  to  write  the 
Hamiltonian  as  the  sum  of  two  terms: 

H • (L  * 2S)  . * I Hk,.  (64) 

K>n,q 
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where  Hkq  is  given  by  equation  (52)  and  the  magnetic  crystal  field 
parameters,  Nkq*  are  given  by  equation  (8).  The  quantity  Hr,  is 
related  to  the  twofold  parameters  N2q. 

Considering  a ferromagnetic  below  its  Curie  point,  we  may  assume 
that  all  the  dipoles  in  the  lattice  are  equal  and  pointing  along  the 
same  direction.  In  particular,  we  assume  a lattice  whose  point  group 
syimietry  is  $4;  thus,  the  only  nonzero  magnetic  crystal  field  parameters 
are  N20,  N4o>  Neo*  ^44*  and  Ne4.  We  assume  the  field  points  along  the 
z-direction.  If  the  material  is  insulating  and  there  are  several  inert 
atoms  between  the  paramagnetic  ions,  we  may  neglect  the  exchange  inter- 
action and  assume  that  there  is  a pure  dipole  interaction.  Under  these 
conditions,  we  may  assume  the  form  of  equation  (64)  for  the  Hamiltonian. 

We  consider  the  zero-field  splitting  of  doublets  by  the  Hamiltonian 
equation  (64).  Usually  only  the  first  term  of  equation  (64)  is  taken 
into  account;  that  is,  it  is  assumed  that  the  splitting  of  doublets  in  a 
ferromagnetic  material  is  due  to  an  effective  magnetic  field  Heff. 
However,  experimentally  it  has  been  shown  that  a s'*ngle  effective  field 
could  not  account  for  all  the  splittings  observed  (say,  by  optical 
methods  or  excited  state  electron  paramagnetic  resonance— EPR) . The 
splittings  can  be  accounted  for  if  all  the  terms  of  equation  (64)  are 
considered.  These  terms  represent  the  effect  of  the  inhomogeneity  of 
the  magnetic  field  near  the  paramagnetic  ion. 

In  conclusion,  we  have  investigated  the  effect  of  a lattice  of 
external  dipoles  on  the  energy  levels  of  a paramagnetic  ion.  Considera- 
tion of  such  an  effect  may  be  represented  by  an  effective  Hamiltonian 
which  is  the  sum  of  two  terms;  the  first  term  representts  the  effects  of 
a constant  field,  and  the  second  term  represents  the  effects  of  the 
inhomogeneity  of  the  field.  No  attempt  has  been  made  to  investigate 
fully  the  consequences  of  the  Hamiltonian  equation  (64)  for  this  report. 


i 
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APPENDIX  A.— CLEBSCH-GORDAN  COEFFICIENTS  AND  RACAH  COEFFICIENTS 
A-1.  CLEBSCH-GORDAN  COEFFICIENTS 

The  Clebsch-Gordan  (C-G)  coefficients  are  defined  through  the 
relation 

* I <ji(mi)j2(m2)|j(ni)>|jimij2m2>  , (A-1) 

which  describes  the  coupling  of  two  angular  momenta  and  ^2  to  ^orm  a 
sum  j.  The  C-G  coefficient. 


<ji(mi)j2{m2)|j("')>  * 


Is  zero  unless  mi  + m2  = m and  unless  |ji  - jal  ^ j 1 iJi  J2I  • som 
ii  + j2  + j must  be  an  Integer. 


Certain  symmetry  equations  relate  different  C-G  coefficients: 

Ji+Ja-j 

<ji(mi)j2('n2)  1 j(n»)>  = (-) 


(A-2) 


= (-)  <ji(-mi)j2(-m2)lj(-"')^ 


(A-3) 


= {-) 


Ji-mi 


2j+l 


<ji(mi)j(-m)  I j2(-n'2)  (A-4) 


= ^ |4^  <j(-m)j2(m2)|ji(-mi)>  . (A-5) 


A-2.'  RACAH  COEFFICIENTS 

In  the  coupling  of  three  angular  momenta,  we  consider  two  coupling 
schemes : 

Scheme  A:  + jj  * J12*  J12  + J3  = J . 

-►  -►  -►  -► 

Jl  + j3  = ji3«  ji3  + j2  * ^ • 
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Scheme  B: 


(A-6) 

{A-7) 
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Coupling  scheme  A Is  represented  by  the  wave  function 

m m m 1 ji2ln’i+m2}><ji2lmi+m2)j3lm3;  | jim)> 

I jlfli  j2'>'2j3''''3>  (A-8) 

and  scheme  B by  the  wave  function 

|B>  = I ^Ji("'i)j3('"3)|Ji3(n'i+">3^><Ji3l'ni+ni3)j2(m2)|j(m)> 

mim2m3 

^ I j2''l2j3''''3^  • (A-9) 

The  coupling  schemes  A and  B are  connected  by  a unitary  transformation 


1B>  = J <A|B>  1A>, 


(A-10) 


and  the  coefficients  of  the  unitary  transformation  are  determined  by 
taking  the  inner  product  of  equation  (A-8)  with  equation  (A-9). 

We  define  the  Racah  coefficients  as  follows: 

W(Ji  jl2  jl3  jl  j)  = [1237"'+  l)(2jj3  + l)f.<A|B>.  (A-11) 

Thus, 


[(2ji2  + l)(2ji3  1)J*  W(j2  ji2  ji3  h'y  Jl  j) 

= I <ji('’’i)j2(n’2)  1 ji2(nii+m2)><ji2(mi+m2)j3(m-mi-m2)  1 j(m): 

m ^(1)2 

<Ji('ni)j3('"-n’i-'n2)  I ji3(m-m2)><ji3(m-m2)j2(m2)  I j(ni)>  . 

The  following  equation  can  be  obtained  from  equation  (A-12): 


(A-12) 
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with  similar  relations  for  the  other  triads.  An  explicit  formula  for  the 
"6j"  symbol  is 


APPENDIX  A 


which  Is  used  in  the  main  body  of  the  report  to  evaluate  explicitly 
certain  "6j''  symbols  and  Racah  coefficients.  In  the  above,  t ranges  over 
all  values  for  which  the  arguments  of  the  factorials  are  nonnegative 
Integers.  The  quantity  &(abc)  Is 


A(abc)  = 


[a  + b - c) ! 


b + c)!  (-a  + b + c)! 


[a  + b + c + 1)! 


(A-18) 


Another  relation  which  Is  used  In  the  text  Is 


The  X-coeffIclent  or  "9j"  symbol,  arising  through  consideration  of 
the  coupling  of  four  angular  momenta.  Is  given  In  terms  of  "6j''  symbols 
by 


A-3.  SOME  USEFUL  RELATIONS 

It  Is  useful  to  represent  the  dot  and  cross  products  of  two  vectors 
In  the  form 

(A  . B)  = ^ (-)^  A^  B_^  . (A-21) 


(A  X B)  = - ^ i I <l(v)l{vi-v)(l(u)>  A^  B^_^  , (A-22) 

1 ^ V 

where 

i Aq  = A^, 

A.,  =^(A^.1Ay),  (A-23) 
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We  list  below  some  useful  relations  involving  the  manipulation  of 
certain  spherical  tensors.  First  of  all,  a third  spherical  tensor  can  be 
formed  from  two  others: 

^ <ki(qi)k2(q2)|k(q)><ki(0)k2(0)|k(0)>  . (A-24) 

kq 

Second,  a relation  for  the  gradient  of  a radial  function  multiplied  by  a 
spherical  tensor  can  be  obtained: 

(grad)^C|^qf(r)  = ^ <l(y) k(q)  | k' (q+y)><l(0)k(0)  1 k' (0)>C,^ . 

[|7  + 4^(r')[k(k+l)(2k+l)]*'^  W(klk'l;  kl)].  (A-25) 

If  f(r)  = r , equation  (A-25)  becomes 

(grad)^Ckqr''={2k+l)<l(y)k(q)|k-l(q+u)><l(0)k(0)|k-l(0)>Ck.i^q+/'''^ 

=-[k(2k+l)fM(u)k(q)lk-l(q+y)>C,^_j^q^^r'^‘^  . (A-26) 

If  f(r)  = r"*^"\  then  equation  (A-25)  becomes 

Q 

(grad)^  f(r)C^q  = -[(k.l)(2k+l)f*<l(M)k(q)  | k+l(q+p)>  . (A-27) 
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APPENDIX  B.  — IRREDUCIBLE  TENSORS 


Irreducible  tensor  operators  are  defined  as  operators  that  have  the 
same  transformation  properties  under  rotation  as  do  the  spherical  har- 
monics. To  express  this  definition  conci‘‘ely,  we  write  the  commutation 
relations  of  the  spherical  tensor  T.  and  the  components  of  the  angular 
momentum 


[J^.  Tkq]  = /knrrrj-  <k(q)  l(w)|k(q  + u)>  , (B-1) 


where  jg  = j^,  Jj  = - yJl  (J^  + i j^),  and  J_j=  (J^  - i Jy) . 

Equation  (B-1)  defines  T^.  as  an  irreducible  tensor  operator  of  rank  k 
and  projection  q. 

The  power  of  the  irreducible  tensor  formalism  lies  in  the  ability  to 
use  the  Wigner-Eckart  (W-E)  theorem  to  calculate  matrix  elements: 


<J'M'lTkq|JM>  = <J(M)  k(q)lJ'(M')><J-llTkl|J>  . (B-2) 

This  important  theorem  states  that  the  dependence  of  the  matrix  element 
on  the  left  of  equation  (B-2)  on  the  projection  q and  the  azimuthal 
Quantum  numbers  M and  M'  is  contained  entirely  in  the  Clebsch-Gordan 
(C-G)  coefficient.  The  so-called  "reduced  matrix  element,"  <J'||Tb||J>. 
depends  only  on  J',  J,  and  k. 

^kiqi  ^kzqz  irreducible  tensors,  then  so  is  the  quantity 

V ' "k2.q-qi  . (B-3) 

This  irreducibility  may  be  shown  by  using  the  W-E  theorem  and  the 
recoupling  schemes  of  appendix  A. 

We  now  proceed  to  calculate  reduced  matrix  elements  of  two  types  of 
irreducible  tensors  in  a basis  |JMLS>:  a tensor  in  L-space  and  a tensor 
in  J-space. 

For  a tensor  in  L-space  (that  is,  a quant;ity  which  satisfies 
equation  (B-1)  with  J replaced  by  t),  we  have,  by  the  W-E  theorem. 
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<J'M'L'S'iT.  (JMLS>  = I <L(M,  )S(M-M,  )|J(M)> 

X <L'  IM|_'  )S'  (M* 

X <L(Mj^)k(q)|L'(ML‘)><L|  IT^I  1L>  , (B-4) 

where  we  have  used  the  expansion 


1JMLS>  = I <L(M|_)S{M-Ml)1J(M)>|LM|_S,  . (B-5) 

If  we  eliminate  the  sum  over  by  the  Kronecker  delta,  recouple  the 
last  two  C-G  coefficients,  and  sum  over  M|_,  we  arrive  at  the  result 

<J'L'S' 1 |T|^|  1JLS>  = [(2L'+1)(2J+1)]^  W(kL'JS'  LJ')  <L'1|T,^||L>  . (B-6) 


For  a tensor  in  J-space,  we  use 
T|(igi  operates  in  L-space,  and  U|^2q2 


the  form  of  equation  (B-3),  where 
operates  in  S-space.  therefore,  we 


<J’M'L'S'|V.  |JMLS>  = I <L(M,  )S(M-M,  )|J(M)> 

^ ^ • 

X <L'(M|^)S'(M'-M|_)|J'(M')><ki(qi)k2(q-qi)|k(q)> 

X <L(ML)ki(qi)|L'(M(_)><S(M-ML)k2(q-qi)|S‘(ML-M(^)> 

X <L'|1T^J|L><S'||U^J|S>  . (B-7) 

If  we  recouple  the  second  and  fifth  C-G  coefficients  and  the  third  and 
fourth  C-G  coefficients  and  sum  over  qi,  we  obtain 
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<J'M'L'S'  IVj^ql  JMLS>=<L'  | | 1L><S'  1 |U^J 

-g-M|  j, 

(-)  ^[(2L'+l)(2S'+l)n2g+l)W(Sk^J'L';S'g) 

X W(LL'kki,;kig)<S(M-M|_)g(M|_+q)|J'(M‘)> 

X <L(-ML)g(ML+q)|k(q)><L(Mj^)S(M-ML)|J(M)>  . (B-8) 

If  we  further  recouple  the  first  and  second  C-G  coefficients  in 
equation  (B-8)  and  sum  over  M.  and  q using  equation  (A-20)  (app  A),  we 
arrive  at  the  final  result;  ^ 

vfL  S 

<J'L'S'llV,^qMJLS>  = [(2L'+l)(2S'+l)(2s+l)(2k+l)]*|ki  kz  k_| 

X <L'||T,^il|L><S'||U^2llS>  . (B-9) 
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